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We present a Ginzburg-Landau theory for a two-band superconductor with empha- 
sis on MgB2. Experiments are proposed which lead to identification of the possible 
scenarios: whether both a- and vr-bands superconduct or a-alone. According to the 
second scenario a microscopic theory of superconducting MgB2 is proposed based 
on the strongly interacting cr-electrons and non-correlated 7r-electrons of boron ions. 
The kinematic and Coulomb interactions of cr-electrons provide the superconduct- 
ing state with an anisotropic gap of s*-wave symmetry. The critical temperature Tc 
has a non-monotonic dependence on the distance r between the centers of gravity 
of a- and 7r-bands. The position of MgB2 on a bell-shaped curve Tc (r) is identified 
in the overdoped region. The derived superconducting density of electronic states is 
in agreement with available experimental and theoretical data. It is argued that the 
effects of pressure are crucial to identify the microscopic origin of superconductivity 
in MgB2. Possibilities for Tc increase are discussed. 



The discovery of superconductivity in MgB2 [1] poses many interesting and fundamental questions 
regarding the nature of the superconducting state as well as the possibility of multi-band super- 
conductivity. The crystal belongs to the space group PQ/mmm or AZi?2-structure where borons 
are packed in honeycomb layers alternating with hexagonal layers of magnesium ions. The ions 
Mg^+ are positioned above the centers of hexagons formed by boron sites and donate their elec- 
trons to the boron planes. The electronic structure is organized by the narrow energy bands with 
near two-fold degenerate u-electrons and the wide-band 7r-electrons. Without any of the lattice 
strain, the a dispersion relations are slightly splitted due to the two boron atoms per unit cell. The 
corresponding portions of the Fermi surface consist of coaxial cylinders along the F - A symmetry 
direction of the Brillouin zone (BZ), whereas the tt- bands are strongly dispersive. In the follow- 
ing sections we present a Ginzburg-Landau (GL) analysis of the two-band superconductor with 
emphasis on MgB2, we apply it in pressure experiments which potentially distinguish the different 
superconducting bands [2], then we provide a microscopic model. 



We distinguish different possibilities: 

Case I: Both a and n bands superconduct. The GL free energy functional for MgB2 in this case 
can be written as : 



(i) If r = then Eq. 1 is the free energy for two bands without pairing transfer (Josephson 
coupling) between them. The quartic term which mixes the two order parameters reflects the 
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fact that although the two bands are different, there is a constraint coming from the common 
chemical potential which connects them in equilibrium. It is also the realization and a measure of 
the strength of the interband interaction of quasiparticles and it affects the critical temperatures. 
The absence of a bilinear in the two order parameters term facilitates the observation of the 
Leggctt's mode in the Josephson tunneling between a two-band superconductor and an ordinary 
superconductor [3] . The onset of the superconducting state in one band does not imply the onset 
in the other. To investigate this possibility in MgB2 an experiment which checks the splitting of 
the transitions due to the lattice deformation by the strain fields can be decisive. In the MgB2 the 
compression due to pressure is anisotropic [4-6]. According to Ref. [4] the compressibility along 
the c-axis is almost twice larger that of the plane compressibility. Therefore application of uniaxial 
pressure on single crystals will affect differently the two gaps and the result will be a splitting of 
the critical temperatures, even if at ambient pressure they appear to be the same. This can be 
measured through a specific heat experiment under pressure. Following Ozaki's formulation [7,8], 
we add to the GL functional Eq.(l) the term which couples the order parameters, in second order, 
with the strain tensor e to first order. 

Fstrain = -Cl{Ti)[d{e^^ + Cj,!,) + e^^] | Vaj^ - C2(ri)(e^^ + Cyy + e^^)|V7r|^ (2) 

where Ci.2(ri) are coupling constants, S is given in terms of the elastic constants of the material 
S = (cii + Ci2 — Ci3)/(ci3 — C33). The effect of the pressure on these two order parameters, each 
one belonging to an one-dimensional representation is to shift the critical temperatures. The new 
critical temperatures can be found, by solving the coupled equations for \ipa,TT\'^ and then setting 
the coefficients of IV'^.ttP to zero. The result is : Tca,v = T^air ~ V<7,-7t x P> where for uniaxial 
pressure 7],^ > ijcr > 0. In order to detect the splitting, the pressure must be above a certain value 
Pmin which dcpeuds on the resolution of the experiment, so that the difference in the two critical 
temperatures AT^ = \T^^ — T"^ + Va- ~ VvlPmin is experimentally detectable. 

(ii) If r ^ then the pair transfering term is present and it means that the onset of superconductivity 

in one band implies automatically the appearance of superconductivity in the other. There is a 
single observed T^ with a different pressure dependence. Analyzing the equations which result 
from the minimization of the GL free energy we get in the regime > T > T^^ : 

aj,aa = (3) 

This gives the pressure dependence of the T^. 

Tcip) = I [T°„ + - + riM + \ - + - + a^] (4) 

and c? = 4r^/(a°a^). Deviations from a straight line at moderate values of pressure can be 

attributed to the two bands. 

Case II: Only the a- bands superconduct. In this case there is one order parameter initially and 
the GL free energy functional is the usual one. The dominant physical situation is an important 
change in the electronic properties of the material under pressure, because the band ((T2) which is 
below the Fermi level at ambient pressure [9] can, partially, overcome the barrier and get above 
the Fermi level at a certain value of pressure (crossover pressure), restoring the degeneracy of the 
two a bands at point T. Then there will be a crossover from a superconducting state ijji to a state 
V'l + V'2- To understand this effect we need to write the GL free energy taking into account the 
second order parameter of the same symmetry as well as the hexagonal symmetry of the boron 
layers. The irreducible representation of the Deh group are four one-dimensional ones (three of 
them have line nodes) and two two-dimensional ones [8] . There are experimental data compatible 
with s-wave or an anisotropic s*-wave order parameter, and a penetration depth experiment [10] 
which suggests nodes or anisotropic .s*-wave. Therefore we consider the gap to be a function of 
and ky with a smooth modulation in the k^ direction as in Ref. [15]. As a consequence the most 
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promising candidates for the initial order parameter are the basis functions of the Fj representation 
of the Dg/i group. To construct the functional of the free energy for two order parameters, wc need 
to take into account the decomposition of the terms containing the derivatives 11, which belong to 
the T5 representation. Then the decomposition 

r5*ori*or5®ri = ri + r2, (5) 

contains the Fi representation once and suggests that there is precisely one term quadratic in Ilx,y 
which mixes the two order parameters (to first order each). There is also a quadratic term which 
mixes the two components, according to the trivial decomposition Fi Fi and a quartic term 
as well. Taking into accounts all the usual terms, we obtain the following expression for the free 
energy: 

F = [ ciM77^|nVl|'+ai|V'l|'+/3l|V'l|' + 77^|nV2r + a2|V2|'+/32|^2|^ + 

71 (V'i>2 + C.C.) + 72(n,^in:v^2* + nj,^in;v2 + C.C.) + /33(|Vi|'IV'2|') + ^^g^^^' }. (6) 

We have used the fact that the 71-term favors the coupling of linear combination of the two 
order parameters with phase difference or tt between them [11], therefore a term V'l*^'! + c.c. 
is incorporated into the /33-term of the free energy. The resemblence between Eq.(l) and Eq.(7) 
is apparent; in the first case the mixing terms are a consequence of the interactions between the 
two bands, in the second case it is due to the modification of the existing gap in the electronic 
spectrum. 

With the present approach we arc in the position to describe the crossover and to consider the 
possibility of different effective masses and parameters (the inclusion of the effect at pc). The effect 
of pressure can be taken into account as usual [12] with the modification due to the particular 
physical situation by writing a2 = a2(T — Tc{p))Q{p — Pc) , 7i = 7i(p^ l)Q(p^Pc)- The additional 
term which couples the order parameter with the strain tensor e, in the regime where p > Pc is : 

Fstrain = -C(ri)[(l + S){€xx + €yy) + 2e^^] X [V'l + ^'2!^. 

The coupling to strain affects both order parameters in the same way. 

The approximate value of the crossover pressure can be estimated as follows. The energy difference 
between the two subbands is approximated by {1 — ricr)'^ A, where the fraction of superconducting 
electrons is I ~ ^ 0.03, i.e. the carrier density per boron atom in Mgi32. The approximate value 
of the crossover pressure Pc which suppresses the deformation potential A can be estimated by the 
expression Pcfl (1 — n<j) \/A where Cl = 30A^ is a unit cell volume of MgS2 and A = 0.04ey^ 
is the deformation potential for a boron displacement u ^ 0.03A [13]. Using these parameters, 
we get a crossover pressure of Pc ~ 30GPa. This estimation shows, that an applied pressure at 
a realistic value influences drastically the electronic structure, producing the degeneracy of two 
bands which arc initially splitted. Superconductivity in the second band will occur at lower values 
of the estimated Pc- The requested band overlapping around the Fermi energy will also occur at 
lower values due to the corrugation of the Fermi surfaces and the already existing strains in the 
material and the anisotropic compressibility. These considerations make the above estimation an 
upper limit of Pc- The degree to which shear stresses of the sample and the surrounding fluid under 
pressure affect the data of pressure measurements is still under investigation [14]. The physics of 
this crossover is similar in spirit to an electronic topological transition of 5/2 kind [16]. In Fig.l in 
a model calculation, we illustrate the expected behavior of T^ and the form of the order parameter 
as a function of pressure at fixed temperature and schematically the kink at Pc of Tc. The anomaly 
close to Pc can be detected directly in a penetration depth experiment under pressure. 
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FIG. 1. The form of Tc as a function of pressure for case II, Pc is the crossover pressure where the second 

band becomes superconducting as well. Inset : The normalized absolute value of the order parameter 
V'(p)/V'(0) = lV'ib)+V'2(p)l/|-!/'i(0)| at T = 0.7Tci. The chosen parameters are: ai = 2(T/Tci-l-0.1p/pc), 
Q2 = (T/Tci - 0.8 - O.lp/pc), /3i = /32 = 1, 133 = Q{p/Pc - 1), 7i = 0.4(p/pc - l)e(p/pc - 1). 



The picture we provide for the second case is complimentary to the first. Both cases may be realized 
since the effects can be detected at different values of pressures. In case I, the value of pressure 
necessary to see the effects is limited cither by the experimental resolution or by the validity of 
the Eq. 4 at low pressure. The specific symmetry of the order parameters is only used to produce 
the additional terms due to strain. In case II there is a larger predicted value for the crossover 
pressure and there are more specific assumptions regarding the symmetry of the order parameter. 
Early pressure experiments [17] demonstrated the overall decrease of Tc with pressure which was 
attributed to the loss of holes. In two of the samples there is a linear dependence of Tc on pressure 
and in two others a weak quadratic dependence. We stress that the samples are polycrystalline and 
the experiment is effectively under hydrostatic pressure. Also the degree of nonstoichiometry was 
not known. The almost linear dependence for a wide range of pressures, makes the GL functional 
as presented, valid for the MgB2. Experiments on single crystals will be able to verify the effects 
which are described. We do not attempt at the moment any fitting of experimental data because 
there is no experimental consensus on the different values of key parameters of the theory (e.g. 
there is a wide range of published data on the value of dTc/dp [14]). 



II. MICROSCOPIC MODEL 



The presented microscopic model of superconductivity in MgB2, is based on the correlated a- 
and the non-correlated 7r-band. The negatively charged boron layers and the positively charged 
magnesium layers provide a lowering of the 7r-band with respect to the u-band which was also 
noticed first in band calcidations [18,19]. In this microscopic treatment the cr-bands arc treated as 
degenarate. The quasi-2D a-electrons are more localized than the 3D 7r-electrons. This leads to an 
enhancement of the on-site electron correlations in the system of degenerated <T-electrons. They 
arc taken to be infinite, while the intersitc Coulomb interactions and clectron-phonon interactions 
simply shift the on-site electron energies. After carrying out a fermion mapping to X-operators 
the Hamiltonian becomes: 
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+ E4"Hp) Kip)^s{p) + H.c] - mE (0 + (0] • (7) 

Here, the X-operators describe the on-site transitions of correlated a-electrons between the one- 
particle ground states (with spin projection s — ±.) and the empty polar electronic states (0) 
of the boron sites and /j, is the chemical potential. The wide 7r-band is shifted with respect to 
the cr-band by an energy r, comprising the mean-field 7r-electron interactions and the electron- 
phonon interactions. We included in Eq. (7) also the nearest neighbour Coulomb repulsion V of 
(T-electrons, which is essential for the low density of hole-carriers in MgB2. The mutual hopping 
between a- and 7r-electrons is assumed to be negligible due to the characteristic space symmetry 
of the orbitals involved. A schematic view of the model is presented in Fig. 2. 
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FIG. 2. A schematic view of the energy band diagram. The degenerate cr-electrons are represented by 
a lower correlated band. 



The diagonal X-operators {X' = X") satisfy the completeness relation, -|- 4X^ = 1, and their 
thermodynamic averages ((A"°), (A"^)) are the Boltzmann populations of the energy levels of the 
unperturbed on-site Hamiltonian in (7). Due to the orbital and spin degeneracy the u-electrons 
occupy their one-particle ground state with density no- = 4 (AT*) per boron site. The correlation 
factor for the degenerated cr-electrons is / = (A"*' -|- A") = 1 — 3 {X'^) = 1 — jUa- It plays an 
important role in all numerics starting from their unperturbed (zeroth order) Green's function, 

D^^ (lu) = // {—iujn — ^J■), whereas for 7r-electrons i?!"'' (w) = 1/ {—iujn ^ 1^ ^ f)- Note that for the 
conventional Hubbard model the correlation factor in the paramagnetic phase is 1 — n/2 (see [23] 
and Refs. therein). 

The energy dispersion of both bands, ^ (p) = ft (p) — /i and e (p) = Sq {p) — r — /i, are gov- 
erned by the zeros of the inverse Green's function {<^,p) = diag {ZJgT^ i'^iP) \ D~^ (w,p)} = 

diag I ; —iojn + £ (p)|, which follows from the Dyson equation i^iP) = (w) -I- 

t {p) , to first order with respect to the tunneling matrix t (p) =:diag{t (p) ; £o (p)}- 
The chemical potential p, for the MgBz =Mg++B- (p2)^=Mg++B- (39"-p"")2-system (Eq. (1)) 
obeys the equation for the total electron density per boron site + 7i;r = 2, with the partial 
electron densities: 

na,n = 2T E e'^^D,,. P) = n.,^ (r, p) . (8) 

The electron densities satisfy the requirements Q < n„ < 1 (due to the correlation factor f, providing 
the quarter-fold narrowing of the degenerate cr-band) and < tt-tt < 2. For any energy difference r 
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between the a- and 7r-bands the chemical potential has to be such that the constraint n^+n^^ = 2 
is satisfied. 

From the system of Eqs. of the constraint and (9), for a flat density of electronic states (DOS 



hereafter) pcr.Tr (&) — (l/2wi,2) ' 



respectively) we derive the chemical potential of the A^^B [p 

5r 



ui2~ with half-bandwidths wi and W2 for a- and 7r-electrons, 

P"'')2 systems as 



W2 



■Wl. 



(9) 



The non-correlated 7r-electrons play the role of a reservoir for the a-electrons. In the energy 
dispersion <^ [p) for the a-electrons, the correlation factor / can be expressed also via the parameters 
Wl, W2 and r as / = {'iwi + W2+ 3r) / (5wi + 4w2). 

The anomalous self-energy for the cr-electrons (Fig. 3) is written self-consistently as 



E(p) = r^ro(p,g). 



s(g) 



e (g) 



V 

s(g) 



(10) 



where the vertex Fq is determined by the amplitudes of the kinematic and Coulomb interactions 
such that Fq (p, q) — — 2tg + V [p — q). We do not include the other kinematic vertices in Fq which 
are essential at a moderate concentration of carriers [23]. 
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FIG. 3. The anomalous self-energy S (p) for cr-electrons. The solid line is an anomalous Green's function. 

In momentum space the Coulomb repulsion between the nearest neighbours (Eq. (7)) reflects the 
tight-binding symmetry of the boron honeycomb lattice [22] . Near the T — A line of the Brillouin 
zone the Coulomb vertex can be factorized as V {p — q) = 2j3t (p) t (q), where the parameter /3 = 
V/6t^ expresses the Coulomb repulsion for the nearest cr-electrons and the energy dispersion is 

t{p) = m(i- 



-y2 — )■ Considering the explicit form of the vertex Fq in Eq. (10), and after 
summation over the Matsubara frequencies ujn — (2n 4- 1) ttT one obtains 



hv)^Yt (.) (1 - ip)) ^ (.) ^^^^yeiq)^^Hq)/2T^ 

'e (q) + S2 (q) 



(11) 



The search for a solution in the form E (p) — TiQ+t (p) Ei converts Eq. (11) for the superconducting 
critical temperature and the gap to 



l^J2*ip)il-f3tip)) 



(12) 



with the gap function E (p) = [1 ~ /3t (p)] Eq. 

Setting the gap function in Eq. (12) equal to zero, one can derive analytically Tc in a logarithmic 
approximation: 
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dWi +4lD2 



+W2-r) (wi + 4r) exp (-i) , 



^ _ (5^»i+4u;2)(3+5/3^;i) / . x 



f3wiW2—2wi —W2 



(13) 



Under the prefactor in the square root the restrictions for an energy shift r guarantee the assumed 
volume of the correlated cr-band, namely Ua- > {r < Wi +W2) and ria- < 1 {r > —wi/A = —t). 
Tc (r) is plotted in Fig. 4 for different values of the parameter (3, reflecting the suppression of 
superconductivity with an increase of the Coulomb repulsion. In the MgB2 case, Tc = 40 K, 
corresponds to r = 0.085 eV and a dimensionless value j3wi — 8 V/3wi. 
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FIG. 4. The non-monotonic dependence of Tc (r) at different magnitudes of the Coulomb repulsion 
(parameter I3w\). Here W2/'W\ — 8. MgB2 is marked in the inset with a solid circle positioned ai r — 0.085 
eV and j3wi = 1.2 for wi — 1 eV and W2 = 8 eV. 



The superconducting gap equation follows from Eq. (12) taken for T=0 

de 



1 = / (e) e (1 - /3e) 



(14) 



It defines the anisotropic superconducting order parameter of s*-wave symmetry. 
The near-cylindrical hole-like cr-Fermi surfaces in MgB2 gives room to calculate the superconduct- 
ing DOS p {E) =J2p^{e- ^e^ (p) + (p)) , where the gap function is S = Eq (1 - /3i {p)) = 

h{l + acos^ i9), where a — (3wi/ (12 (1 — (3wi)), 6 = Sq (1 ~ P^i) and ^ is the azimuthal angle. 
Then the DOS normalized with respect to the normal DOS is 



E 



dz 



E^ - {1 + az'^y 



(15) 



For a Coulomb repulsion such that (3wi < 1 the parameter satisfies a > and the superconducting 
DOS becomes 



p{b < E < {l + a)b) 
piE>{l + a) b) 



5f,T(sin-iy^ 



(16) 



which is expressed in terms of the complete and incomplete elliptic integrals K and T, respectively 
with modulus q = \/ (E ~ b)/2E. 
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FIG. 5. The superconducting DOS, normalized with respect to the normal DOS, for the Coulomb 
parameter range /Jwi > 1. At /Swi = 0.92, the result by Haas and Maki [15] is reproduced (inset). 

The DOS (16) has cusps at E=zL (1 + a)b. A similar result was obtained in Ref. [15] for a non- 
specified parameter a > 0. In our case the parameter a is controlled by the in-plane Coulomb 
repulsion V as the authors of Ref. [20] noted. At l3wi = 0.92 the DOS of Eq. (16) (see Fig. 5, 
inset) reproduces Fig. 1(b) of Ref. [15]. 

For our case of an enhanced Coulomb repulsion f3wi > 1, we have to take the parameter a < in 
the gap function and the superconducting DOS (Fig. 4) is then given by 



p((l-|o|) 6 <£:<&) 
Po {E) 

p{E>b) 
PO (E) 



y/iE+b)\a\b 



2\a\b-^ 



E+b 
2E 



(17) 



It contains two logarithmic divergencies at E — ±b and a gap in the energy range \E\ < (1 — |a|) b. 
The two-gap ratio is 1/ (1 — |a|). 

Measurements on MgB2 with scanning tunneling spectroscopy [24,25] and with high-resolution 
photo-emission spectroscopy [26] revealed the presence of these two gap sizes. From the ratio 3.3 
between the two gaps in Ref. [24] we can extract the parameters |a| « 2/3 and /3wi « 1.14, whereas 
from data of Ref. [25] one can derive Pwi = 1.21. A value (3wi = 1.14 can be estimated from 
Ref. [26]. Point-contact spectroscopy [27] shows gaps at 2.8 and 7 meV, from which we estimate 
a Coulomb repulsion parameter f3wi — 1.16. The recent study of energy gaps in superconducting 
MgB2 by specific-heat measurements revealed two gaps at 2.0 meV and 7.3 meV [28] for which 
(3wi = 1.13, and a gap ratio 3 — 2.2 [29], for which f3wi = 1.14 — 1.15. Measurements of the 
specific heat of Mg^^B2 also give evidence for a second energy gap [30]. Raman measurements [31] 
established pronounced peaks, corresponding with gaps at 100 cm~^ and 44 cm~^. From these data 
one can extract a = —0.56 and Pwi = 1.18. At Pwi = 1 we have gapless like superconductivity. 
In this case the superconducting DOS is linear with respect to the energy near the nodes of the 
superconducting order parameter. Then the superconducting specific heat Ce ~ and the NMR 
boron relaxation rate ~ at low temperatures. From this point of view it is interesting that the 
data of Ref. [32] shows a Ce ^ behaviour and a deviation from the exponential BCS behaviour 
in Tf ^ (T) of [33] visuahzcd in MgBa. 



III. CONCLUSIONS 



We have done a GL analysis of the superconductivity of a two-band superconductor with par- 
ticular attention to MgB2. In that framework, we provide the different possibilities for a T-P 
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phase diagram, make predictions for the role of the bands and discuss different experiments from 
which crucial information can be extracted. Then we have analyzed the superconducting prop- 
erties of the material MgB2 within the framework of a correlated model Eq. (7). The existing 
electron-phonon and non-phonon approaches to the superconducting mechanism in MgB2 can be 
separated in two groups: one pays attention to the cr-electrons and the other to the 7r-clcctron 
subsystem. We have taken into account both the correlated a- and noncorrelated 7r-electrons. 
Analysis of our results leads to the conclusion that superconductivity occurs in the subsystem of 
(T-electrons with degenerate narrow energy bands whereas the wide-band 7r-clcctrons play the role 
of a reservoir. Superconductivity is driven by a non-phonon kinematic interaction in the cr-band. 
A lot of evidences in favour of two different superconducting gaps can be explained by anisotropic 
superconductivity with an order parameter of s*-wave symmetry, indiiced by the in-plane Coulomb 
repulsion. For an enhanced interboron Coulomb repulsion {j3wi > 1) the logarithmic divergencies 
in the superconducting DOS (Eq. 17, Fig. 5) are manifested by a second gap in the experiments. 
In our approach the electron-phonon coupling is hidden in the parameter r . Therefore the pressure 
and isotope effects can be explained by the dependence of r (w) on the phonon modes. From the 
non-monotonic Tc dependence on r it follows that the MgB2 material is in the underdoped regime 
(around r > —wi/A). For a fictitious system A^+B2; where the two electrons arc contributed by 
atom A, the superconducting critical temperature increases with an r increase. A pressure increase 
lowers the a-band with an r decrease resulting in a negative pressure derivative of Tc in agreement 
with experiment [34]. The bell shaped curve Tc{r), with the MgB2 position in the underdoped 
regime, shows a possibility to reach higher T^'s in diboride materials with an AIB2 crystal struc- 
ture. We suggest the synthesis of materials with increased r- values (e.g. with "negative chemical 
pressure") and optimized smaller interatomic B-B distances in the honeycomb plane. 
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